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Based on the metastable electron-pair energy band in a two-dimensional (2D) periodic potential 
obtained previously by Hai and Castelano [J. Phys.: Condens. Matter 26, 115502 (2014)], we present 
in this work a Hamiltonian of many electrons consisting of single electrons and electron pairs in the 
2D system. The electron-pair states are metastable of energies higher than those of the single¬ 
electron states at low electron density. We assume two different scenarios for the single-electron 
band. When it is considered as the lowest conduction band of a crystal, we compare the obtained 
Hamiltonian with the phenomenological model Hamiltonian of a boson-fermion mixture proposed 
by Friedberg and Lee [Phys. Rev. B 40, 6745 (1989)]. Single-electron-electron-pair and electron- 
pair-electron-pair interaction terms appear in our Hamiltonian and the interaction potentials can be 
determined from the electron-electron Coulomb interactions. When we consider the single-electron 
band as the highest valence band of a crystal, we show that holes in this valence band are important 
for stabilization of the electron-pair states in the system. 
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The mechanism of electron pairing in unconventional 
superconducting materials still remains as one of the ma¬ 
jor and unresolved problems in condensed matter physics. 
It is usually believed that the 2D characteristic (such 
as a layered crystal structure of the cuprate supercon¬ 
ductors) and strong electron-electron correlation are im¬ 
portant for unconventional superconductivity. From 
the experimental observation of the very small coherence 
length in hight T c superconductors Friedberg and Leejjl) 
(referred to as FL below) proposed a phenomenological 
theory, the so-called s-channel theory. In this theory of 
the boson-fermion model, they introduced phenomeno¬ 
logically a local boson field (j) representing a pair state 
which is localized in coordinate space. Each individual 
(j) quantum is assumed to be unstable with an excita¬ 
tion energy 2v. It has been shown that this assump¬ 
tion makes it possible for the s-channel theory to exhibit 
many superconducting characteristics0- f _9]. This model 
has been studied both for three-dimensional (3D) and 
two-dimensional (2D) systems. As a matter of fact, the 
boson-fermion model has been extensively studied in the 
area of unconventional superconductivity. [Ml It can 
be traced back before the BCS theory for conventional 
superconductivity. M 

The mechanism of electron pairing and the nature of 
the bosonic excitations are still key issues for the under¬ 
standing of unconventional superconductivity. [3] Based 
on previously obtained metastable electron-pair states 
in a 2D crystal[l5|, we present in this work a Hamilto¬ 
nian of a 2D electron system consisting of single electrons 
(fermions) and electron pairs (bosons). This new Hamil¬ 
tonian obtained from electronic structure calculations of 
a 2D crystal includes single-electron states, electron-pair 
states, as well as Coulomb interactions between single 
electrons and electron pairs. The interaction potentials 
present in the Hamiltonian can be obtained from the elec¬ 
tronic structure calculations with electron Coulomb in¬ 


teraction. 

We start with the single-electron states in a 2D square 
lattice crystal. The crystal periodic potential is given by 
V c (x,y ) = Vo[cos(g*) + cos(<?y)], where q = 2tt/X with A 
the period and Vo the amplitude of the crystal potential. 
The Schrodinger equation for a single electron is given 
by H(r)?/i„ i k(i") = -E^kVvk(r) with the single-electron 
Hamiltonian H( r) = —V 2 + V c (x, y) in the units of effec¬ 
tive Bohr radius a b and effective Rydberg R y , where k is 
the wavevector in the first Brillouin zone, n is the band 
index, E n k and tpn, k(r) are the eigenvalue and eigenfunc¬ 
tion, respectively. 

When we consider two interacting electrons in this pe¬ 
riodic potential, the Hamiltonian is given by H( iq, r 2 ) = 
H(ri) + IF(r2)-|-2/|ri — r 2 |. The ground state of this two- 
electron system can be found for |ri — r 2 | —» oo due to 
the electron-electron repulsion. It means that the ground 
state consists of two non-interacting single electrons sep¬ 
arated by an infinitely long distance. However, it was 
found previously that metastable electron-pair states [jj| 
may exist in this system. In such an electron-pair state, 
two electrons of spin singlet stay in the same unit cell of 
the crystal lattice in the relative coordinates r = rq — r 2 
with an average distance between them only about one 
third of the lattice period. Because two electrons in the 
same unit cell are of maximum overlap, they are strongly 
correlated. Strong electron-electron correlation reduces 
greatly the Coulomb repulsion between them. When the 
amplitude Vq of the crystal potential is larger than a 
certain value, strong correlation between two electrons 
together with the crystal potential modulation leads to a 
metastable electron-pair state. The electron-pair states 
in the center of mass coordinates R=(ri+r 2 )/2 are given 
by Bloch wavefunctions and form energy bands. The 
eigenenergies of the electron pairs are given by E^ au and 
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the wavefunctions can be written as 
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where G; = l x qi + l y qj (with l x ,l y = 0, ±1,±2,...) is 
the reciprocal-lattice vector, R nm (r) (with n=0,l,2,..., 
and m= 0 , 2 ,..., n) is taken from the eigenfunction of a 
hydrogen atom with a scaling factor of Laguerre function 
as given in Ref.|l5j, bo = 2 and b m = 1 for m > 1. It was 
found that the lowest eigenvalue of an electron pair is spin 
singlet and forms an energy band. In terms of the energy 
per electron, this electron-pair band is found in between 
the two lowest single-electron bands. The minimum of 
the electron-pair band at the T point in the first Brillouin 
zone is higher than the maximum of the lowest single¬ 
electron band (at the M point) in the studied square 
lattice potential in Ref. jl5l | . Their difference is defined 
as an energy gap E g . 


In the presence of many electrons in the crystal, we 
expect that for a certain electron density many-body in¬ 
teractions may stabilize the electron-pair states leading 
to a macroscopic electron-pair density. Without loss of 
generality, we consider two energy bands: a lower band 
E | of single-electron states V'k(r) and a higher band 
v o + E^ air /2 (given by energy per electron) of electron- 
pair states. As was shown in Ref. [Hi], the minima of the 
two bands in the square lattice are at the T point in the 
first Brillouin zone and their difference at low-electron 
density is given by vq = E g + Ebw > where Ebw is the 
band width of the single-electron band. 


In order to obtain the many-electron Hamiltonian, we 
will consider two different scenarios for the single-electron 
band. We consider first the single-electron band as the 
lowest conduction band in a real crystal. Free conduction 
electrons occupy this band. It is empty in an undoped 
insulator at low temperature. When the electron-pair 
band is taken into account, some electrons may form 
metastable pairs in this band. We assume that there 
are N electrons consisting of Np single electrons and 
Nb electron pairs in the system (N = Np + 2Np)- At 
low electron density, all conduction electrons occupy the 
single-electron band, i.e., N = Np and Nb = 0. At 
higher densities, electron pairs may appear and many- 
body interaction may stabilize these pairs. The Hamil¬ 


tonian of the system can be written as, 


JV N N 

H= Y,H 0 (r l ) + -J2J2 J — 

i=l 


2 zz ' iTj — Tj | 

i= 1 j^i ' 1 


Nf 


1 Nf N f 9 


i=i £ i=i&i' T ' 


Nb 

E 

ot—1 


E ff °( r “> 1 ') + oEE 


i 


U=1 

Nb Nb 2 2 

EEEE 


2 ^ / 1 ^OL.V 1 * 0 ;, a\ 

1/=1 1 ’ 1 


a—1 (3ty=L n 
Np Nb 2 




i —1 ct=1 u=l 


r, - r n 


(1) 


where r.^ denotes single electron i , and r a ^ electron pair 
a (with two electrons i/=l and 2). The first two terms in 
the above Hamiltonian are Np single electrons and the 
Coulomb interactions between single electrons. The next 
two terms are Nb electron-pairs and the Coulomb inter¬ 
actions between the pairs. The last term is the Coulomb 
interactions between single electrons and electron pairs. 
Due to possible transitions: two single electrons an 
electron pair, Np and Nb are not constants. Two single 
electrons can form an electron pair through the Coulomb 
interaction between them. In the opposite process, an 
electron pair can decouple into two single electrons. A 
macroscopic number of electron pairs in the system de¬ 
pends on the electronic band structure of the crystal (de¬ 
termined by the crystal potential), the total electron den¬ 
sity and temperature, etc.. When there are no electron 
pairs in the system, the above Hamiltonian reduces to 
that of interacting single electrons with N = Np. 

Using the eigenfunctions of the single-electron and 
electron-pair states, we can obtain the field operators of 
single electrons and electron pairs. The single electron 
(fermion) field operators can be written as, 

vU r ) = E^ k W Ck > ff 

k 


and 


^( r ) = E^k(r)ciL> ( 3 ) 

k 


where the operators cj^. and Ck. CT are creation and an¬ 
nihilation operators, respectively, for a single electron 
of momentum Ek and spin a. They obey the fermion 
anti-commutation relations {ck,cr, c^, a ,} = ^k.k'^o-.cr', 

{Ck,<T,Ck',<7'} = 0, and {cl a ,cl, a ,} = 0. 

The electron-pair (boson) field operators are given by 

T(R,r) = ^T k (R,r)6 k (4) 

k 
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and 


*t(R )r ) = 5>*(R ir )6t, (5) 

k 

where the operators b £ and b k are creation and anni¬ 
hilation operators, respectively, for an electron pair of 
momentum ftk. They obey the commutation relations 
[bk, &£,] = 4,k', [6k, 6k-] = 0, and [6^, b\.,] = 0. 

Considering that the electron-pair number Nb = {N — 
Nf)/2 in the system depends on many-particle interac¬ 
tions and band re-normalization, the Hamiltonian of the 
system in second quantization with single-electron and 
electron-pair states can be written in two parts, 


H = H 0 + H 1 . 


( 6 ) 


The first part Hq corresponds formally to the five terms 
in Eq. ©, given by 


Ho = 


+ 

+ 


F -6'k C k.(T Ck : fT 


k,<r 

l 

2 A 


E E V q C k 1 -q.a C k 2 +q,a' C k2,<?' c k 1 ,<T 


ki,k2,q o-.cr' 


E [ B £“ r + H ^ 


E ”pj>(9) i, k 1 -, i> k,+„ ll lo i, ki 
ki,k 2 .q 

YX E vS ep^) b lx- q C k+ q ,Ai C k,a, ( 7 ) 

k,ki,q,cr 


where v q =2— is the electron-electron Coulomb poten¬ 
tial, Vpp(q) = 4 v q fp P (q) the pair-pair interaction poten¬ 
tial with a form factor f pp (q), and vl p (q) = 2 v q f ep (q) the 
single-electron-electron-pair interaction (scattering) po¬ 
tential with form factor f ep (q)- Notice that, in Eq.® we 
have taken the bottom of the single-electron band as the 
reference for energy. The summation over q does not in¬ 
clude q=0 because it cancelled out with the background 
ion-ion interaction and the system is neutral. The form 
factors fep(q) and f pp {q) reflect finite size of the electron- 
pair (against point charge) and are given by 





q (<) 


FIG. 1: The form factors for the electron-pair-electron-pair 
(fp P (q)) and single-electron-electron-pair ( fe P (q )) interaction 
potentials, and the interaction potential v t (q) for single- 
electrons-electron-pair transition. A = 1.5 an and Vo = 12, 
15, and 18 R y . 
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(9) where J m (x) is the Bessel function of the first kind. 
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The second part Hi of the total Hamiltonian Eq. © is 
responsible for the transition or transformation between 
two single electrons and an electron pair. The interaction 
potential is just the Coulomb potential between the two 













4 


electrons. It is given by 
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where v t (q) is the single-electron-electron-pair interac¬ 
tion potential for transition (transformation) which is 
given by 
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drR nm (r)J m (qr). 


( 12 ) 

The form factors f PP (q) and f ep (q ) and the potential 
v t (q) can be obtained from numerical calculation. They 
are given in Fig. 1 for the square lattice with A = 1.5 as 
and Vo = 12, 15, and 18 R y . 

The total electron number operator is given by 


N = 


r)ijj<j(r) + 2 dR dr¥ (R, r)$(R, r) 


— E c k.CT C k ,<t + 2 E 


(13) 
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It commutates with the total Hamiltonian Eq. © and is 
conserved. 

In comparison with the FL Hamiltonianjj, the interac¬ 
tion terms due to pair-pair and single-electron-electron- 
pair interactions appear in our Hamiltonian given by 
Eqs. ©, (|T|). and (JTTJ). The FL model Hamiltonian of 
the boson-fermion mixture was given by Eqs.(1.7), (1.8) 
and (1.9) in Ref. |§i containing three terms only, which 
correspond to the single electrons (or holes), the single 
bosons (i.e., electron pairs) and the interaction between 
electrons (or holes) and the bosons. The boson quantum 
was introduced phenomenologically with spin zero and 
with twice the mass and charge of a single electron, and 
it was further assumed unstable with an excitation energy 
2u. But, the detailed microscopic structure of the boson 
quantum and its origin of the pairing mechanism were un¬ 
known. Therefore, the interaction potential between the 
electrons and bosons was described through a coupling 
constant as a parameter. In this work, however, we build 
a more complete theory with the Hamiltonian given in 
Eqs. ©, ©, and (fill) based on the electronic structure 
calculations of a two-dimensional square lattice. From 
the spin-singlet metastable states of the electron pairs, 
together with the single-electron states in the system, we 
have obtained quantitatively the interaction potentials, 
such as pair-pair interaction and single-electron-electron- 
pair interaction (both for scattering and transformation), 
presented in the Hamiltonian starting from the electron- 
electron Coulomb potential. 


We consider now another possible scenario for the 
single-electron band. In this scenario, the single-electron 
band is considered as the highest valence band of an in¬ 
sulator. The conduction carriers in this band are holes 
with positive charge and dispersion relation E^. For an 
intrinsic insulator, this band is fully occupied by electrons 
at low temperature and consequently the hole density is 
zero. If we assume that there are Np^ holes (from dop¬ 
ing initially) in the valence band and Nb electron pairs 
in the electron-pair band (and > 2Nb ), the total 

carrier charge number (i.e., the net hole number) in the 
system is given by Nh = N^ — 2 Nb- Transition of two 
electrons from the valence band to the electron-pair band 
creates an electron pair leaving two holes in the valence 
band. The scenario is shown schematically in Fig. 2. 
Taking the top of the valence band as the reference for 
energy, the many-particle Hamiltonian in this scenario 
can be written as, 


H (v > = H, 
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where the operators d £ and dk.o- are creation and anni¬ 
hilation operators, respectively, for a hole of momentum 
fik and spin a. They obey the fermion anti-commutation 
relations {dk.cr; d^, — dk.k‘'^er,er / 5 {^k,<r 5 ^k'} — 0 ; 

and { d l a ,dl, a ,} = 0. The hole-hole interaction poten¬ 
tial v q is the same as that of the electron-electron in¬ 
teraction. The pair-pair interaction potential is also the 
same as that in Eq. 0 , but the hole-electron-pair inter¬ 
action (scattering) potential is different from the corre¬ 
sponding single-electron-electron-pair potential by a sign, 
v h P (q) = ~ v e P (q )> because the hole has positive charge. 
The corresponding form factors are also given by Eqs. ( 8 ) 
and (9). The second part H^ ' 1 of the total Hamiltonian 
is equivalent to H\ given in Eq. (HID . It is responsible for 
creating (destructing) an electron pair in the electron- 
pair band and two holes in the valence band at the same 
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time. From the definition of holes, hole creation (annihi¬ 
lation) is actually electron annihilation (creation) in the 
valence band. So the creation of an electron pair in the 
electron-pair band and two holes in the valence band is 
equivalent to the transition of two single electrons from 
the valence band to the electron-pair band, and the in¬ 
teraction potential v l (q) keeps the same form giving by 
Eq.(12). 

The total (net hole) number operator is given by 

^ = E4a-2E^- ( 17 ) 

k,cr k 

It commutates with the total Hamiltonian Eq. ED and 
is conserved. 



FIG. 2: A schematic plot of the energy band structure with 
electron pairs and holes. 

In this scenario, the stabilization of the electron pairs 
can be understood as following. At low density, the at¬ 
traction between electron pairs and holes are dominant 
leading to the formation of exciton-type states. [H] The 
most possible scenario should be an electron pair binds 
two holes into a biexciton (X 2 ). In two dimensions, the 
energy of an exciton (for m e = nih) is approximately 
E * = — 2R y . The energy of a biexciton is given by 
E= 2 E* + AE^ 2 3 , where A E^ 2 is the binding en¬ 
ergy of the biexciton in relation to two free excitons. In 
a 2D system, this binding energy is about 20% of the 
exciton energy E*. 0 Therefore, the energy of an elec¬ 
tron pair can be reduced with about 4.4 R a . In other 
words, the gap energy E g can be reduced by 2.2 R y . At 
high density, the system tends to be an electron-pair and 
hole plasma as soon as the electron pairs can stabilized 


through many-particle interactions. Because the electron 
pairs and holes cannot bind together due to strong inter¬ 
action or screening of many particles, in this case, the 
band-gap renormalization due to pair-pair and hole-hole 
interactions are essential to stabilize the electron pairs. 
Band-gap renormalization (BGR) has been extensively 
studied in nonlinear optics of semiconductors, where one 
deals with an electron-hole plasma (EHP). In an EHP, 
the BGR is given by a sum of electron and hole self¬ 
energies and is a function of the particle density and 
temperature. jl8l-l20| It is known that the BGR at finite 
density can be several times the exciton binding energy 
in a 2D system. For a high density EHP, the BGR is 
mainly induced by the Coulomb repulsion among the par¬ 
ticles (the Coulomb hole effect). Therefore, we believe 
that the BGR in our present electron-pair-hole system 
should be of similar value, i.e., the renormalization of the 
electron-pair energy will be enhanced at high density and 
can be several times the biexciton energy. Consequently, 
the band-gap normalization in the present system can be 
larger than the single particle gap E g leading to a neg¬ 
ative gap. As soon as the bottom of the electron-pair 
band touches the Fermi energy Ep of the valence band, 
a macroscopic electron-pair density is stabilized in the 
system. 

In summary, we present a Hamiltonian for a 2D elec¬ 
tron system with single electrons and electron pairs. The 
electron-pair states are metastable at low electron density 
and they are expected to be stabilized through many- 
particle interactions. We have considered two different 
scenarios for the single-electron band. In the first case, 
the single-electron band is considered as the lowest con¬ 
duction band. In comparison with the FL Hamiltonian, 
extra terms such as single-electron-electron-pair and pair- 
pair interactions appear in our Hamiltonian. The inter¬ 
action potentials in the obtained Hamiltonian are deter¬ 
mined from the electro-electron Coulomb interaction. In 
the second scenario, the single-electron band is consid¬ 
ered as the highest valence band. We make a connec¬ 
tion of the electron pairing and its stabilization to hole 
density in the system. Our present study has shown ba¬ 
sic insights into possible electron paring mechanism in a 
crystal. 
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